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CN . Abstract 

>■ 

^^ ' In this work we discuss an analytic approach towards the solution of 

OO 
,_^ I pure Yang- Mills theory in 3 + 1 dimensional spacetime which strongly 

^ I suggests that the recent strategy already applied to pure Yang-Mills 

^ ' theory in 2 + 1 can be extended to 3 -|- 1 dimensions. We show that 

^D ■ the local gauge invariant variables introduced by Bars gives a natural 

,i^ . generalisation to any dimension of the formalism of Karabali and Nair 

' , which recently led to a new understanding of the physics of QCD in 

D , dimension 2 + 1. After discussing the kinematics of these variables, we 

'T^ I compute the jacobian between the Yang-Mills and Bars variables and 

.5^ ' propose a regularization procedure which preserves a generalisation of 

p\ ' holomorphic invar iance. We discuss the construction of the QCD hamil- 

j^ ' tonian properly regularized and compute the behavior of the vacuum 

wave functional both at weak and strong coupling. We argue that this 

formalism allows the developpement of a strong coupling expansion in 

the continuum by computing the first local eigenstate of the kinetic part 

of Yang-Mills hamiltonian. 

1 Introduction 

In a remarkable series of papers Karabali and Nair P] and Karabali, Kim 
and Nair j2] (see [3] for an introductory review) have developed a novel and 
powerful method in order to adress the confinement problem of pure Yang-Mills 
in dimension 2+1 (I will refer to this work as KKN for short). Their approach, 
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inspired by the work of Feynman ^ on the Schrodinger representation of QCD, 
has led to an analytic determination of the string tension and the possibility 
to devise a strong coupling expansion in the continuum and evidence for a 
mass gap. The main idea behind this work is to formulate the Yang-Mills 
hamiltonian in terms of local gauge invariant variables. This parametrization 
heavily uses the power of the complex structure which is available in 2D. 

This work has been followed by by a deeper investigation of the vacuum 
wave functional in the Large N limit by Leigh, Minic and Yelnikov |2l IH] . Bas- 
ing their analysis on a quadratic ansatz for the wave functional, these authors 
have investigated the structure of the Schrodinger equation and the action of 
the kinetic term on a restricted class of operators. From this study they have 
proposed an interpolating kernel capturing the behavior of the vacuum wave 
function both in the infrared and the ultraviolet. This kernel, expressed sim- 
ply as a ratio of Bessel functions, exhibits an infinite set of resonances arguing 
in favor of a constituent picture for glueballs which becomes apparent in the 
KKN formalism. This have led them to a theoretical prediction of the glue- 
ball spectra simply related to the zeros of Bessel function which is in striking 
agreement with the large N lattice data [7]. 

This series of work open many new questions and avenues in the study 
of QCD, first in 2+1 dimensions where among others things one can wonder 
about the possibility of exploring along these lines the mesons spectra, the 
inclusion of a non zero temperature and where a deeper understanding of the 
renormalization group in this context is needed. But one of the most pressing 
issue is to know wether such a methodology can be extended and applied to the 
case of pure Yang-Mills in dimension 3+1. The main message of the present 
work is to show that this is indeed the case. 

The first obstacle to overcome lies in the fact that the KKN formalism uses 
the full strength of complex analysis which is available in 2 space dimensions 
but not in 3D. The key and simple insight obtained in collaboration with R. 
Leigh and D. Minic |0 is the fact that the KKN gauge invariant variables 
are related to lines integral from infinity to a point x and admit a natural 
generalization in any dimension which was first proposed by I. Bars j9| some 
time ago and baptized 'corner variables'. 

In this work we explore the consequences of this insight and develop a 
new set of techniques allowing us to work in terms of the Bars variables and 
extract a first set of physical information from this formalism. In section |21 



we explore the kinematics of Bars variables and reestablish some important 
formulae relating them to the usual Yang-Mills variables, we also identify the 
key notion of an 'holomorphic symmetry'. In section |21 we present one of our 
first main result which is the computation of the jacobian of the transformation 
between Yang-Mills and Bars variables. This computation is possible once we 
devise a natural regularization scheme which preserves holomorphic symmetry. 
The jacobian is then shown to be trivial in the case of real coordinates and 
given in terms of a generalization of the WZW action in the case of semi- 
complex coordinates. This result implies that the trivial wave functional is 
normalizable and hence the theory is amenable to a strong coupling expansion. 
This strongly resonates with the results obtained in lattice QCD but in the 
context of a continuum formulation. 

In section |3] we present the construction of the regularized hamiltonian in 
Bars variables and as a check give the construction of the wave functional in 
the ultraviolet regime. 

In section we show that the regularized kinetic term acts diagonally on 
the potential term which is our second main result. In contrast with 2+1 D 
the eigenvalue, having a dimension of a mass, is cutoff dependent. This allows 
us to construct the vacuum wave function at first order in a strong coupling 
expansion and using an argument of dimensional reduction we show how this 
determines the value of the string tension in this regime. We conclude on the 
open issues, one of the main one being the problem of the continuum limit. 

In a joint companion paper |H] we explore the physical consequences of the 
results obtained here and strongly argue in favor of a deep similarity between 
the 2+1 and 3+1 dimensional case. 



2 The kinematics of Bars Corner variables 

We denote by Ai, i = 1,2,- --D the spatial components of the Yang-Mills 
connection. These are taken to be anti-hermitian and our main interest is D = 
3, but we keep it free in the first part of the paper. One of the key ingredient of 
Karabali Kim and Nair (KKN for short) formalism is a parametrization of the 
configuration space of 2d Yang-Mills in terms of local gauge invariant variables 
which uses heavily at first sight the available complex structure of 2 dimension. 
In fact jHI the KKN gauge invariant variables can be understood as Wilson 



lines integral from infinity to a point x and admit a natural generalization in 
any dimension which was proposed by I. Bars [Hj- We denote them Mi{x), 
they satisfy the defining equation 

A = -diM,M-\ (1) 

No summation over repeated indices is assumed here and after. A solution is 
given by 

Mj{x) = exp[- / A] (2) 

where the integral is a straight spatial contour for fixed x* for i y^ j and exp 
is the path ordered exponential, explicitly 

oo „ 

Mi(x) =l + ^(-l)" / A,{ti,X2,X3)---Ai{tn,X2,X3)dti---dtn. (3) 

Gauge transformations 

A ^ Af = gAg-' + gd^g-' (4) 

act on M's as 

Mi ^ gMi. (5) 

One can then define gauge invariant local variables 

H,, = M-'M,. (6) 

Note that Hjj = 1 and Hji = H~-^ -this just means traversing the corner in 
the opposite direction is precisely the inverse element. The Hij are unitary in 
a real coordinate basis. There is also a constraint (here written for D = ?>) 

HijHji = 1, HijHjkHki = 1. (7) 

The translation to 2+1 KN variables is obtained if one uses a complex basis 

M, = M, Ms = {M^y\ H = Hs, = M^M (8) 

Note that if we had used a real coordinate basis, then we would have had 
a pair of (unrelated) unitary matrices Mi, M2. In 3+1 there is no complex 
structure and we will work in a real basis where Mj are unitary matrices. One 



Huz = 


= H 


Hzu = 


= H^ 


Hsz = 


= H^H 



could also choose to work in a "semi-complex" coordinate basis {-u, z, z} for 
D = 3, then we could parameterize this as 

(9) 
(10) 

(11) 

For example, one could use the notation M^ = M, Mz = (M^)~^, MlMu = 1, 
with H = MIM. The constraint takes the form H^uHuzHzz = 1- In other 
words, there is in D = 3 in the semi-complex coordinate basis a complex H- 
field (compared to a Hermitian field in Z) = 2); thus there are twice as many 
degrees of freedom. 

A very important feature of this parametrization is a generalization to any 
dimension of holomorphic invariance, and even if there is no complex structure 
in dimension 3 we will refer to this symmetry as holomorphic symmetry. This 
extra symmetry acts as 

Mi ^ Mihi\x'), J ^ i (12) 

The condition j ^ i on the function hi is the analogue of holomorphy. This 
leaves the gauge fields invariant, and one finds 

H,j ^ KH.jhT^ (13) 

In the complex basis, we would have M^ — > Muh~^{z^z), M^ -^ h{u,z)M\ 
M -^ Mh^{u,z) and so H ^ hu{z,z)Hh^{u,z), W -^ h{u,z)H^h-\z,z). 
Note that one can introduce currents 

J,, = id,H,,)Hr^\ (14) 

(in 2-1-1, J ~ Jzz and J''" ~ —Jzz)- The Jij transform as ('holomorphic') 
connections 

Jij ^ hiJijh^^ + djhi h:[^ (15) 

(This is not a typo - it only depends on hi.) 

In the real coordinate basis, it appears that there are six currents that are 
apparently distinct. However, there is a 'reality' condition on their derivatives 
of the form 

d,J,, = -H,^{d,J,i)H7^ (16) 
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(in D = 2, this reads dJ = H{dJ'^)H ^) By defining J^- = —HijJjiH^j^, we 
may rewrite this as 

and so there are covariant derivatives Dij = dj — Jij. These currents are related 
to the magnetic field Fij = diAj — djAi + [Ai, Aj] by 

d^Ji, = -Mr^F.jMi. (18) 

Note that the covariant derivatives can be written in terms of the usual deriva- 
tive and Mi as 

V, = d, + A = MAM-\ (19) 

Let us finally note that Wilson loops observables can be expressed entirely 
in terms of the currents: Lets consider C a closed curve in the 12 plane. The 
Wilson loop observable in the representation R can be expressed as 

Wr{C) = Tm'^-^c^^dx, = Tr^V^c-^i2(^(t))dx2_ ^20) 

This is easily obtained once we write the connection as a gauge transformation 
of the current 

Aidxi + ^2^X2 = -dxi^iMiMf ^ - dx2(92M2M2"^ (21) 

= -{dxidiMiM^^ + dx2iMiJi2M^^ - d2MiM{^)){22) 

If one works out what this correspondence imply for a rectangular plaquette C 
in the plane 1/3 = cste with corners (0, 0)(i/i, 0)(0, y2)(l/i, I/2) oriented clockwise 
one simply get 

= Tr^ {Hu{yi,0)H^2\yuy2)Hi2{0,y2)H{2\0,0)) (23) 

In order to continue and work out some formulae necessary for the computa- 
tion of the hamiltonian in these variables lets introduce some notations and 
conventions. The connection is expanded in terms of anti-hermitian gener- 
ators Ta (in order to avoid unnecessary factors of i), Ai = A'^Ta satisfying 
the algebra [T)j,T(,] = fab'^Tc, we also denote by Tr the trace in the vectorial 
representation, so that Tr(l) = A^. by — 2Tr(TaTb) = 6ab (the minus sign is 
because of the antihermiticity and the 2 is the standard convention). The 
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Yang- Mills action is taken to be Sym = 2^ / Tr(F^j,F^,^) and the hamiltonian 
in the Yang-Mills variables is 



V^ / 9' f ^W 1 ..a^2 



2, a 



^-i.y-T^M«j +v<^°)^ <2^) 



with Ff = \eijkF°-^. This should be supplemented by the gauss law constraint 

In the adjoint representation the generators are given by (T°-)bc = —fabc, 
where index are raised or lowered with the metric 6ab and the trace in the 
adjoint is denoted by Tra^. A group element M is represented in the ad- 
joint by Mab = -2TT{TaMThM-^), clearly we have {M-'^)ab = M^a and also 
MTbM-^ = TaM^b- 

We denote Gi{x,y) the inverse of di, d^Gi{x,y) = S{x,y) (no summation 
on repeated indices). Since di admits zero modes, its inverse is not uniquely 
defined and we will work with the explicit choice Gi{x) = 6{xi)6{x2)S{x3) 
where 9 is the heaviside function and Gi{x,y) = Gi{x — y) (which is not 
antisymmetric). This choice is not arbitrary, it is the unique choice consistent 
with the definition of the variables Mi as an ordered exponential. Indeed we 
can rewrite Q in terms of the propagator as 



M,{x) =Y,i-^T j MG^A,r{x,y), 



(25) 



where {GiAif{x,y) = J dzGi{x,z)Ai{z)Gi{z,y)Ai{y) etc... 

One can first compute the derivative of Mj with respect to Ai. Starting 
from 

AI = 2Tr (T^d^MiM-^), 

one obtains the relation 

SA'} = -{M,Y^d,{M-^5Mif 

with the obvious notation M^^5Mi = {M~^5Mi)'^Tb. Inverting this relation 
gives 

^ \j^a/l = -S}G,iy,^)iM-\x))\ = 5l{M,G^{x,y)J>. (26) 



where we introduced the notation 

G\x,y)^-G{y,x). (27) 

t stands for transpose and the minus sign insures that d^G*{x, y) = S{x, y). In 
other words we can express Ai derivatives in terms of right derivative on the 
group 

[P:ix),M^iy)]^6lM,iy)TJix,y) 

as 

^^ = {M,{x)\' j dyG\{x,y)Pl{y) ^ {M,G\P\{x). (28) 

The next step is to express A derivatives in terms of the currents, from their 
definitions one has 

Jij = djHijHr.^ = -Mr^ {A, + djMiMr^) M, = -Mr^A^Mi - Mr^djMi, 

taking its variation one gets 

5Jij = -Mr^6A,Mi - [Aj, Mr'6Mi] (29) 

with Dij = dj — Jij. If one takes (j^Uj) into account this reads 

6 J,, = -MrHAjM, + [A,, G^MrHA^Ml (30) 

from which we get 



5A1{x) 



-{M,{x))a'S{x,y), (31) 

{Dl)\G.{y,x){M-\x))\ (32) 



6Al{x) 

= 0, i^j,k (33) 



sJUy) 



6Af(x) 



where 



(note that (Tj^b = -fcab)- 

So far everything is similar to 2D. Before going on, it is useful to dwell 
further on the kinematical structure and write down more precisely the struc- 
ture of phase space in this new variables. In order to do so we need some 



notations: one introduces the momenta variables which replace the electric 
field generators in our variables and the corresponding operator smeared with 
a Lie algebra valued 1-form field ^i{x) = $"(x)Ta 

Ulix)^iGlP%ix), Ui^)^Jdx^ti^)Ki^)- (34) 

Obviously, 

n(a$) = -P($), n($) = -P{G(!>), (35) 

with {d(f))i = di(f)i, (G$)j = Gi^i and P($) the smeared version of P^. The 
algebra in terms of the P, M variables is simply 

[P($), M,{x)]^ = MM^), [P($), Pm, = Pm n- (36) 

Where the bracket index q denotes quantum commutators to be distinguished 
from the Lie algebra commutators. It will be useful to know also the commu- 
tator with the currents 

[P($), J,,]^ = H,,d,cP,H7.' - A,0.. (37) 

The generators of gauge symmetries can be naturally expressed in terms of 
these, if X{x) = X"'{x)Ta labels the parameter of the infinitesimal gauge 
transformation one can write 

Gx = -[tt (xV,^) = P{M,XM-^). (38) 

The generator of the holomorphic symmetry can be written in a similar form 
in terms of Lie algebra elements H^ = H^Ta satisfying diHi = 0, 

Gh = P{H,). (39) 

One sees that the conditions of holomorphic invariant Gh ~ is necessary in 
order for the momentum LI to be well defined since if Hi is an holomorphic 
transformation Il{dH) = = P{H). Moreover the observables LI, J are gauge 
invariant 

[Gx,Pm, = [Gx,J^J], = 0. (40) 

When written in terms these observables the algebra reads 

[n($),n(vi/)]^ = -u{d[G^,G^]), (41) 

[Ui^),Mx)]^ = -{H,,^,H7^){x) + {D,,G.,^i){x), (42) 

[■h^.Jki], = 0. (43) 
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where [^,% = [$„ ^i] and obviously 5[G$, G^] = [$,G^] + [G$,^]. It 
might be interesting to note that this algebra possesses strong similarity with 
a centrally extended algebra, this is clear if one takes the derivative of (j42j) 
leading to 

[n(<l>), d^ij] = [{G^h ^^J^J] + A,$. - di{H,,^,Hr^). (44) 



3 Determinant and regularization 

The wave functions in the Schrodinger representation of pure Yang-Mills are 
gauge invariant functionals of Ai and the scalar product is given by 



11^11'= / Dii{A)^{A)<^{A), (45) 

JA/g 

where the integral is over the space of gauge connections modulo gauge trans- 
formations and Dfi{A) = y^T^r = voucf' -^^^^ ^^^ previous section we know 
that we can equivalently describe gauge invariant wave functionals as holomor- 
phic invariant wave functionals of Hij or Jij. In order to express the physical 
scalar product in this new variables we need to compute the jacobian. Since 
Ai = —diMiMf^, 6Ai = —{yi5Mi)M^^, the change of variables involve a 
determinant 

e^ = det ( ^^' \ = det(ViV2V3). (46) 

The variational derivative of the action is given by 

^^ = Tr.4(V.)-^(x,a;)r.] (47) 

the trace being in the adjoint representation. Since Vj = MidiM~^ the co- 
variant propagator can be expressed in terms of the standard propagator as a 
coincident limit 

-^ = hm {M,{x)Gi{x,y)M7\y))'' f,,^. (48) 

A similar coefficient arises when one expresses the kinetic term in terms of the 
group variables using pHjl 

(49) 
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Integration over x and summations over indices should be understood. The 
relation between these coefficients is not a coincidence, indeed the presence of 
e^ in the integration measure insures that (|49p is a self adjoint operator. 

Of course, in order to make sense of these statements, one needs a regu- 
larization of the propagator which preserves all the symmetries, namely the 
holomorphic symmetry Mj i— >■ Mih~^. Since Vj = MidiM^^ is invariant un- 
der holomorphic transformations this means that the regularized propagator 
should transform as Gf (x, y) ^— i> hi(x)G^{x, y)hl^(y), where fx denotes the mo- 
mentum cut-off scale (this is obviously true for the unregulated propagator 
which satisfies hi{x)Gi{x, y)h~^ (y) = Gi{x,y)). 

The key ingredient in the construction of the regularized propagator are 
the group valued functionals Aj(x, y) defined to be 

Ai{x,y) = Hi2iyi,X2,X3)H23{yi,y2,X3)Hsi{yi,y2,y3), (50) 

and cyclic permutation for A2,A3. These group elements transform under 
holomorphic transformation Hij \-^ hiHijhJ^ as 

Ai(x, y) ^ hi{x2, a;3)Ai(x, y)h^^{y2, 2/3) (51) 

and similarly for A2, A3. They also satisfy the key properties 

dfAi{x,y) = 0, and Ai{x,x) = 1. (52) 

We can now propose a regularization scheme which preserves holomorphic 
invariance : 

Gi{x,y)= dzGi{x,z)Ai{z,y)6f,{z,y) (53) 



where the regularised delta function is given by (say) 

3 
S^.{x, y) = Yl ^m(^» - yi)^ ^m(^») = -^e"'^'^?. (54) 



This regularised propagator satisfies 

afGf(x,i/) = A,(x,y)(5^(x,y). (55) 

Moreover the integral ()53|1 can be explicitly performed and one gets 
Gi (x, y) = Ai(a;, y)6^{xi - yi)5^,{x2 - 2/2)'^/.(x3 - ys) = Ai(a;, i/)G'^(x, y), (56) 
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where 9f^{x) = -k^ J^'^ dte~* is one half plus half the error function. 

Now that we have a regularization which preserves gauge invariance and 
holomorphic symmetry we can evaluate the propagator at coincident points 

2 

and compute Gj'(a;, x) = ^. This propagator is proportional to the identity 
hence ri- = and 

I det(ViV2V3)| = |det(ai 92^3)1 (57) 

the determinant is independent of the connection. This is independent on 
the form of the regulated delta function. We have proven this results here 
in dimension 3 but the formalism is valid in any dimension therefore in any 
dimension D 

|det(Vi---Vz5)| = |det(9i---9D)|. (58) 

Since this is for us an important conclusion, let us give an other totally 
independent proof of the same result using lattice gauge theory. Putting gauge 
theory on a lattice provides a gauge invariant regularization of the theory. We 
start with a gauge theory on a square periodic D-dimensional lattice, we choose 
one origin and an orientation and denote the displacement vector of one lattice 
unit in the direction i by Cj, the sites of the lattice are labeled by a; = ^^ XiCi 
with x^Nti being identified with x. The gauge connection is encoded in terms 
of group elements gi{x) associated with the link {x,x + e^). The analogous of 
the Bars variables are given by holonomies starting at one corner of the lattice. 
Let X = 'Y^^XiCi, < Xi < N, then 

Mi{x) = gi{xY^gl^{x - e^) ■ --gi^ix - XiCi), (59) 

so that the starting element is always on the plane Xi = 0. We introduce the 
discrete derivative dif{x) = /(x) — /(x — Cj), the Mj satisfy the difference 
equation 

diMi{x)M-\x) = (1 - g,{x)), Mi{x) = 1, forx^ = 0. (60) 

which is clearly the discrete analog of (Q). It can be equivalently written as 
VjMj = 0, where the discrete covariant derivative is VjM(a;) = gi{x)M{x) — 
M{x — Cj). Now, we want to compute the Jacobian of the transformation 
gi — > Mi, since gi{x) = Mi{x — ei)Mf^{x) we first remark that Mj depends only 
on gi so that the transformation matrices is block diagonal and the problem 
is essentially one dimensional. A direct computation gives the variation 

giix)-'6giix) = iVi6M,{x))Mr\x) = Mi{x){diMr'6Mi{x))Mi{x). (61) 

12 



Therefore 

D D 



det(Vi---VD)| = \{\dei{Vi)\ = \{\dei{MidiM-^)\ 

i=l i=l 

D 

= J]|det(a,)| = |det(ai---aB)| (62) 



i=l 



Where we have used the block diagonal form of the transformation matrix, and 
also the factorization and invariance under conjugation of the determinants 
which are all valid operations since we are in a finite dimensional context. 

There is an important subtlety here which is hidden in the domain of 
definition of the operator di. This comes from the fact that Mi{x) is not a 
periodic function on the lattice, since on one side of the lattice Mj(0) = 1 and 
on the other side Mj(iVej) is arbitrary. This means that the determinant of 
di which arises from the change of variables should be computed in the space 
of Lie algebra valued functions satisfying (^(0) = and (f){Nei) arbitrary and 
not on the space of periodic Lie algebra valued functions. One can easily see 
that then the residual derivative determinant can be exactly evaluated and is 
in fact equal to unity. This simple answer can be traced back to the lattice 
regularization which starts from group variables and doesn't affect our main 
statement which is that the full determinant is independent of the connection. 

3.1 Gauge invariant measure 

This results allows us, up to a gauge field independent determinant, to make 
the change of variables from Ai to Mj. The induced measure VMi on the fields 
Mi should be thought as being the product over spacetime points of the Haar 
measure DM = Yl^. dM{x). 

More precisely, the proper way to define the measure is to give its moments 
evaluated on a proper class of functionals. The functionals which have a 
well defined integration under this measure are called cylindrical functions 
F{M{x)) = F„(M(xi), ..., M(x„)) which depends on the value of the field 
only in a finite number of points and the measure DM is defined by its value 
on cylindrical functions to be J DMF{M{x)) = J dMi...dMnFn{Mi, ..., Mn) 
the integral being over a product over Haar measure. The Haar measure is 
right-left invariant and one can write (since the measure is ultralocal, in the 
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sense just defined, we can look at only a fix x) 

dMidM2dM3f{Mi,M2,M3) = [ dMidHi2dHisf{Mu MiH^, M^His) 



(the order of integration being important). If / is gauge invariant the integral 
over Ml factorises, supposing the measure on the unitary group to be nor- 
malized we get J dHi2dHi3f{Hi2,Hi3). This shows that the measure on A/G 
is 

JDHuDHrs = J DH,2DH,sDH2s5,{H,2H,sH23) 

where in the last equality we insert a delta function on the group to empha- 
size the symmetric form of the measure under permutation of indices. The 
interpretation of this delta function constraint should be clear: it is the in- 
tegrated version of the Bianchi identity e^^^ViFjk = expressed in terms of 
gauge invariant observables. 

Now, the identity functional ip{H) = 1 can be viewed as a limit of identity 
cylindrical functional the integral of which being always one if one choose the 
normalized Haar measure. So we can conclude from this analysis that the 
trivial wave functional is integrable. 

This result is purely kinematical but it captures an essential feature of the 
formalism, namely the fact that even if we are working in the continuum, the 
choice of Bars variables allows us to recover results which are easily obtained in 
the lattice formulation, and as we will see, allows us to devise a strong coupling 
expansion. If one computes the expectation value of Wilson lines observables 
in this trivial kinematical vacuum one sees, for rectangular Wilson loops (j23j) 
and from the definition of the measure, that 

{WR{y)) = j DH,2DH,3WR{y) = 6r,o- (63) 

Wilson loop expectation values are non zero only for trivial representation, 
this confirms the interpretation that this trivial vacuum corresponds to an 
infinite string tension, which is what one expects in a crude g —>■ oo limit. 
Now, one could say that such a result seems at odds with the well known fact 
that in 2D there is a non trivial determinant in the KKN formalism, namely 
the hermitian WZW action mi21- This apparent contradiction can be easily 
resolved since, thanks to the work of Gawedzki and Kupianen ^Uj, we know 
that there is a deep relationship between hermitian and unitary model. 
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More precisely if we specify their results to the genus zero case, the relation 
between correlation functions of WZW theory level k with SU{N) group and 
field g{x) and the hyperbolic model with field h{x) = m(x)"''m(x) is given by 






4=1 n=l 

where A/y is the number of conformal blocks (independent of Xi and the sphere 
metric). When extrapolated to fc = clearly this identity shows the normal- 
izability of the trivial wave functional. 

Now it is also well known that the correlation functions of the unitary WZW 
model are zero unless the highest weight K^ is integrable, that is (A/j, 9) < k 
where 6 is the highest weight of the adjoint representation {6, 6) = 2. If 
we extrapolate these results to A; = this means that only primary fields 
associated to the trivial representation and its descendant are non zero in 
agreement with (jM|) . Moreover as noticed by KKN |2] the correlators of the 
currents of the unitary WZW model can be obtained by analytic continuation 
{k + A^) —>■ —{k + N) from the hermitian ones. The reason being that both 
satisfy the KZ equations with opposite {k + N) and both are invariant combi- 
nations under monodromies, which fix them uniquely up to normalization. 

3.2 Semi-complex coordinates 

The determinant is trivial in the real basis where H are unitary matrices. 
However, if one choose to work in the semi complex basis we find a non trivial 
determinant which is a generalization of WZW model in higher dimension and 
is related to the magnetic mass term introduced in ^T]. For completeness we 
present this computation here. 

We have seen in the first section that in the semi-complex coordinates 
X = (x, X, Ux) the gauge invariant data is entirely encoded into a complex 
SL{N, C) field H = Huz- The unregulated propagators are given by 

G„(X,y) = e{ux-Uy)5\x-y), (64) 

Gz{X,Y) = ^ M ux-Uy), (65) 

7r(x — z) 

G-ziX.Y) = —^-—6{ux-uy). (66) 

nix — z) 
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The integrals (fS^ defining the regulated propagators can be explicitly per- 
formed leading to 

GC(X,r) = 9,iu,-Uy)5l{x-y)AuiX,Y), (67) 

G.iX,Y) = ^'I'l^Z'j^ (l - e-^^l-^I^A.(X,y)) , (68) 

G,(X,F) = M!f^^(l-e-^'l--^l%(X,r)), (69) 

Ti[x — y) y / 

with |a;p = xx and 

A„(X,r) = H{x,x,Uy){H^Hr\y,x,Uy)H^{y,y,Uy), (70) 

A,(X,r) = (/7t/7)-i(^,a;,^jijt(^^^^^^)^(^^^^„^)^ (71) 

A,{X,Y) = H\x,y,u,)H{x,y,Uy){H^H)-\y,y,Uy). 
The value of the propagator at coincident point is given by 

G(:(^,^) = g, (72) 

G,(X,X) = -MH^H)-'d{WH), G,{X,X) = Md{H^H)){WH)-\ 

The variational derivative of the determinant (j48j) can be explicitly computed 
and after some algebra one gets 



^r = -^ fTr{{H^H)-^6{H^H)d{{H^H)-'dH^H)} 

112 J 



(73) 



This equation can be easily integrated out in terms of a three dimensional 
generalization of a Wess-Zumino-Witten action 



3 , 

7r2 _ 1 



r = i fdud'^zTi {d{H^H)d{H^H)-') + ^ fduf Tr {[{H^ H)-^d{H^H)]^) 



2Nfi 
The last integral being over a three ball Bu bounding the plane u = cste. 

4 Hamiltonian 

In this section we write down the Hamiltonian in terms of gauge invariant 
variables. From now on and in the following one shall stick to three (space) 
dimensions, a real basis and we introduce some notations well adapted to this 
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case. The vacuum wave functional we look for, is expressed as a functional 
of the currents J12, J23, Jsi- These should be thought as a vector (or three 
dimensional holomorphic connection) and we denote 

Hi = i/j-ij, Ji = Ji-ii = diHiH^ , Di = Di^ii = di — Ji (74) 

and we introduce the 'magnetic field'^ 

5i_i = diJi+i. (75) 

(all indices are modulo three indeed) The potential term is readily expressed 
in these variables 

V = -Y, j^^{{^^J^+lf{x))dx = -Y, j'^<B,{x)fdx (76) 

i i 

i,a 

and with the help of the holomorphic regularization one can compute the ki- 
netic term ()49|1 : the coefficient of the linear term vanishes for reasons identical 
to the vanishing of the variation of the jacobian and we are left with 

r = -^ E l{GlP:y{x){GlP:){x)dx = \Y.j PMQfiy, ^)Pt{^)dydz 

i i 

(78) 
Here for later convenience we have introduced a (unregularized at this stage) 
kernel 



Qf{y,z) = 5'^\G,G])iy,z) = S'^" J G,iy,x)G,iz,x)dx. (79) 

and we have used that (G-P^)"''(x) = {P^Gi){x). 

If one look at the detailed structure of the kinetic term one encounters 
a disturbing infrared divergence hidden in it, namely if one look at ()78|) the 
structure of each term is given by 

-- / d^xdyidzie{yi,Xi)e{zi,Xi)P^{yi, X2, Xs)P^{zi,X2, X3) (80) 



^It is related to the true magnetic field Fi by i?i-i = — Af^ ^Fi^iMi 
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1^ 



using the identity 1 = 9{x,y) + 9{y, x) we can express this term as 

-]^j dxAxJj "Pfdyi') (j "pfdZiVx2,X3) 

+ - j d^xdyidzie{yi,xi)e{xi,zi)P'^{yi,X2,x^)P^{zi,X2,x^). (82) 
The first term in this expansion contains 



Pt{yi,X2,x,)dy,= / Tr(Pi(y)C(y))dy 

which is a generator of an infinitesimal holomorphic symmetry with generator 
0"(l/) = —2T"'5{x2 — y2)S{x3 — ys), di(j)'^{y) = and therefore acts trivially on 
holomorphic invariant states. We can therefore drop this term and consider 
the kinetic term to be given by the second term ()82j) . This amounts to replace 
in all the previous expressions G)j = GG* by the convolution product 9, = GG. 
This second term is now free of any infrared problem, the integral over xi can 
be performed and we are left with 

-- dyidzidx2dx3\yi- zi\TT{Pi{yi,X2,X3)Pi{zi,X2,X3)) (83) 

which exhibits a stringy nature of the kinetic term with some linearly rising po- 
tential. There is even a deeper justification arguing also in favor of this second 
form of the kinetic term which comes from the derivation, in the connection 
variables, of its matrix elements given by 



We also want our theory to be CPT invariant which means that instead of 
conjugating with simple complex conjugation \E' ^ ^ we can equivalently re- 
place this conjugation by a CPT transformation ^ —>■ CPT{^). The CPT 
transformation of the gauge field is given CPT{Ai{x)) = —A}^{—x) = Ai{—x). 
When we change variables from Ai to Mj we need to integrate over one di- 
mension from — oo to x, this breaks parity and the CPT conjugate of Mj is 
not just the hermitian conjugate : it is given by a path ordered integral which 
starts from +oo that is CPT{Mi{x)) = {Mi{-x)hiy. Where hi{x2,X3) = 
■^J_oo i{-xi-x2-x3)dxi ^ |g ^j^g holonomy across space at fixed X2,xs, and dihi = 
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0. Of course if the wave function is holomorphic invariant it doesn't depend on 
hi and CPT invariance is satisfied. Now if one starts from a form of the scalar 
product with is manifestly CPT invariant (that is CPT{^) is used instead 
of ^) and compute from there the form of the kinetic term one can see that 
the A derivative acting on the right can be expressed in terms of the M and 
the propagator G but the derivatives acting on the left should be expressed 
in terms of the propagator G* and in this case the form of the kinetic term is 
expressed in terms of = GG. This form of the kinetic term being the one in 
which CPT invariance is manifest. That is the kinetic term we will use is 

^ = ^ E / n{ym\y. z)Pl{z)AyAz (85) 

i,a,b 

Here the unregularised kernel is 

Qt{y,z) = 5''\G.^,){y,z) = S'^' f G,{y,x)G,{x,z)dx. (86) 

This kernel satisfies 

d^Qfiy, z) = 5-'G,{y, z), d^Qf{y, z) = -5-'>G,{y, z), (87) 

d^d-tQf{y,z) = -5-H{y,z). (88) 

We can now express the kinetic term in the current variables, first recalling 
for comfort and according to our new notations the commutator 

[Pa{^\ JKv)] = SHH,{y)fadl6{x,y) - 6f^\Dl^,)\6{x,y). (89) 
Therefore 

r=i^/d,d.n-te..)^^ (90) 

with 

^t{v.^) = [{D^f''{D^rQ^-i{y,z)-6'%y,z)], (91) 

^r^uiy,^) = \[Dl,Hr\z)f{G,{y,z)-G,{z,y)). (92) 



The full QCD hamiltonian is 



n = g^T + ^V. (93) 

9 
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4.1 Regularised Hamiltonian 

In order to define mathematically the theory we need to include a regulator. 
Preserving holomorphic invariance is a strong constraint on the form of the 
regulator and we have seen that there is a natural holomorphic invariant reg- 
ularization with scale parameter /i which is available in our context. This 
regularization amounts to replace 6°'''Gi{x, y) by its holomorphic invariant reg- 
ulated version (jKH|) Gf " (x, y) , the holomorphic regularization of the kernel 
is denoted Q'^ and is given by 



r^{x,z) 



et{y,z) ^ {G^G^ny,z) = J dx {Uy,x)Ux,z))"' G^iy,x)G';, 

(94) 
with 

Giiy, x) = e^ijji - xi)5^{y2 - X2)5^{y^ - xg), 

and (see section IH)), 

Ai(y,x) = H2{xi,y2,y'i)H^{xi,X2,y'i)Hi{xi,X2,x:i) (95) 

plus cyclic permutations. Under holomorphic transformations we have 

Mi^Mih-\ H,^hi_iHih~\ P^^ (hi) a' Pi,, Bi_, ^ h^Bi.ihr' 

Ai{y,x) ^ hi{y)Ai{y,x)hi{xy^ (96) 

and Q^iy,z) ^ h,{y)e>;{y,z)K{z)-\ 

The regularised kinetic term written in a manifestly hermitian form is given 
by 

^ = i E /(Gf POi(x)(GfP^(x) + {G^ni{x){GfP%{x)dx, (97) 

i,a 

= \J2 J Pi^(y)(^f^(y^') + ^"^'('^y))Pc{z)dydz, (98) 

i,b,c 

i -^ 
where we used the transposed kernels 

{Gfr\x,y) ^ -(Gn''"(2/,x), {efr\x,y) ^ {eff'^{y,x). 
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It is convenient for our subsequent computations to move all the derivatives 
to the right; in doing so we encounter a commutator 

[P^{y),Q^'%y, z)] ^ hm [PUy), Q^'^y, z)l (100) 



fi—*oo 



where Pi is a left invariant derivative regulated according to our prescription: 
it satisfies (GfP*)^ = {GiP^)a and is explicitly given by 

Pl^abj) = ldzAf{y,z)S^{y,z)Pl^{z). (101) 

The calculation of this commutator gives 

[PUy),etiy,^)]dy = -<dri^,z), (102) 



where G^ = O^^/ab^- This term is of order 0(-^). For the other commutator 
we find 

j[n{y).Qt%,y)]dy = p'i{zl (103) 

with 

p^{z) = ldx{G>:{z,x)G^{x,z,^)r (104) 

where z^-^ = (xi, 2:2, -Zs), z^^ = {zi,X2,Z3) etc... This term is also of order O(^) 
and independent of Zi, it corresponds to an holomorphic transformation. We 
can now write the kinetic term in terms of the currents : in order to do so one 
needs to evaluate the double commutator 

et{y,z)[Pl,{y),[P:{z),r,{x)]]dydz (105) 



which is also of order O ( ^ ) . The conclusion is that up to terms / dzil'ii 



'6Jf{z) 



with VL1{z) = O i\\ the regularised kinetic term takes the form (plj) with 

n'^{y,z) = [DyD^Q,^,{y,z) + H,{yMdtQ,{y,z)){H,{z))-'f, (106) 
^luiy,^) = -l[Dl,Hr\z)dt {e.iy,z) + el{y,z))]'\ (107) 

where {e^)''\x,y) = e''"(z,x) and {DDQf^ = D\D^dQ'"^. 
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5 Vacuum wave functional 

In this section we show how one can, in our formahsm, determine the behavior 
of the vacuum wave function both in the weak and strong couphng regime. 

5.1 Weak coupling regime 

Let us first look at the weak couphng hmit of the theory : as is usual one 
first rescale the currents J -^ J/ g and then expand the hamiltonian in power 
series in g,H = ^^ g^Ti-i. After field rescaling, the potential term V is of order 
zero, therefore this expansion can be realized entirely as an expansion in the 
kinetic term and more precisely in the coefficients Q. This expansion is an 
expansion in the number of currents J. Since J has a dimension of mass it 
always enters the expansion in the combination J/p, if one think in terms of 
Fourier modes. This means that we should interpret this expansion as being 
valid in the ultraviolet sector p — » oo of the theory and Tio is the Hamiltonian 
governing the gluon dynamic in the deep ultraviolet, this Hamiltonian being 
free. 

The diagonal terms Qa contains only a finite number of terms up to g"^ 
but the off diagonal component due to the presence of the group element H 
contains an infinite number of terms in this expansion. We are now interested 
only in the structure of the first non trivial term Tig and the coefficients of the 
kinetic term at zero order are given by 

fiP(l/,^) = [{^ym)e^^^{y,z)-6{y,z)]6'^, (108) 

^Tuiy,^) = l[^UG^{y,^)-G,{z,y))]6'^. (109) 

This expression is not particularly illuminating at first sight, however if one 
writes down this hamiltonian in Fourier space B{p) = J dxe~^P^B{x) and 
takes as fundamental variables the magnetic fields Bi = (9j+i Jj_i, after simple 
algebra the hamiltonian becomes 



I f d^p 
rin 



'^J C^ 
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'-(2.f i/S, - ,.,,) j4) J^ + Sf (rtBr(-p) 



The Poincare invariance of the theory becomes manifest. The matrix con- 
tracting the derivatives has eigenvalues p^,p^,0 the eigenstates associated to 
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p^ correspond to the two gluons polarization perpendicular to p. The eigen- 
state corresponding to the eigenvalue is puzzling at first sight since we are 
working in a gauge invariant formalism. In order to properly understand it one 
needs to remember that the gauge invariant measure constructed previously 
contains only two scalar independent variables per Lie algebra generator, this 
is encoded in the constraint HiH2H^ = 1 or in differential form in the Bianchi 
identity Vj-F* = 0. This means that Ji, J2, J3 are not independent variables, 
and that any operator proportional to this constraint is strongly 0. When 
expressed in terms of our local variables and in the limit g ^ this constraint 
imposed by the measure reads J2iPi^iip) — 0- 

We can now easily diagonalize the free hamiltonian Tio subject to this 
constraint and give explicitly the ground state in this regime 

nB) - exp -— J ^^ (110) 

which is indeed the one expected and where we have reinstalled the dependence 
in g by rescaling to the original variables. 

5.2 Evidence for confinement 

The previous result is not very impressive by itself since it just shows that 
our formalism can handle appropriately the perturbative regime which is well 
understood by other means. However, if one look back at the formulation in 
terms of Bars variables this formalism, even if it is local, is very similar to 
the lattice formulation of QCD which is easily amenable to strong coupling 
expansion. It shows that the route we have chosen allows us to deal also with 
the strong coupling or infrared regime of the theory. 

In the strong coupling regime g ^ 00 the hamiltonian H is dominated by 
its kinetic term, as emphasized for instance in the hamiltonian formulation of 
lattice QCD [T2|. At first order in the strong coupling expansion the vacuum 
state is just the constant wave functional which is normalizable as we have seen 
even if we work in the continuum. In order to go beyond this naive result and 
construct the vacuum state in the infrared regime, we establish as a first key 
result that the potential term V viewed as a wave functional is an eigenstate 
of the kinetic term. Since the eigenvalue has a dimension of a mass we see the 
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emergence of a dynamically generated mass scale, namely 

/T- : V := M : V : with M = ^, (111) 

2(27r)i 

where N is the dual coxeter number N6^ = f^'^'^ftpq (with our convention on 
Lie algebra generators) the :: denotes normal ordering of the operator and ■ 
denotes the action of the operator T and /i is our ultraviolet cutoff regulator. 
It is important to keep in mind that this result is valid only in the strong 
coupling regime, the true QCD mass scale should appear in the continuum 
limit as the regulator is removed /z — > oo. 

The first step of the computation follows directly from the definition of the 
regulated Hamiltonian (j^ I106|) which is defined in terms of the holomorphi- 
cally invariant regulated Of (y, -z)"^ and its derivative. From this definition, 
one obtained after integration by part and noticing that the off diagonal terms 
do not play any role 

g^T- f -TT{d,^iJi)\x)dx (112) 
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,bb 



dy. 



y=z 



The RHS of pi2|l is the sum of two terms which can be computed separately, 
this relatively lengthy calculation is presented in detail in the appendix. The 
computation of the first term of ()112|) gives 

- [CiSf-iTr.. (BfBfe^ife.B))]^.^ = (113) 

jJi^[-2iVT.(Bi.)-V(iV=-l)]+0(i). 

The computation of the second term of ()112j) gives the same result ! 

- [CiCiTra. {HM{dyd!ef{y,z)){H,{z)r')]^^^ = (114) 

/^ r riArrp/7-)2 N A 2/Ar2 iNl , /^ I ^ 



[-2NTriBl,)-Af^'iN'-l)]+o(- 



4(27r) 

The factor N in front of the trace is the dual coxeter number coming from 
the relation between trace in the adjoint and in the vectorial TTad{T°'T^) = 
2NTt{T"-T''). Summing over i shows that ()112|1 is 
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Therefore if one defines the normal ordered operator 

/V2 _ 1 

:V:=V-6/i2^^, (116) 

then one gets the announced result pilj) . It is interesting to note that T acts 
diagonally not only on the Lorentz invariant combination V but also on each 
term in the sum that is 

N^ -1 
r- : Tt{B^) := M : Tr(5f ) :, : Tr(S,^) := Tr(S,f) + 2/i^^^. (117) 

This result is extremely simple and strongly resonates with the 2 dimensional 
results. It suggests that the action of the kinetic term on local gauge invariant 
operators is analogous to a kind of dilation operator which counts the number 
of J constituents of the operators it acts on. In order to validate such a picture 
higher order computation should be performed. 

5.3 strong coupling expansion 

We can now exploit the full strength of our formalism and show that the 
strong coupling expansion can be performed following a strategy devised in 
2+1 by KN and strongly reminiscent of the cluster expansion of Lattice QCD 
[T!^ . The QCD hamiltonian reads (lets take a finite dimensional analog for 
the presentation of the argument) 

H = g^T+ 4t V, with T = n'^did^ (118) 

9 

Then we see that in the strong coupling regime g —>■ oo the hamiltonian reduces 
to its kinetic term and the trivial wave function \1' = 1 is a solution in this 
regime to the equation Hip = up to order l/g"^, this solution is normalisable 
as we have seen in the previous section. In order to go further lets suppose 
that 

n^^5iSjV = MV, (119) 

which is what we have proven in the last section in the context of QCD, and 
lets consider the wave function \1/ = e 9^ $, using the previous equality we 
have 

He~^'^^ = e"^^^$, (120) 
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where 

1 y 

and we see that the wave function $ = 1 or ^ = e a^ is a vacuum solution 
at order 1/ g'^ and therefore describes the vacuum in the strong couphng regime. 
The higher order terms in this expansion can be recursively constructed in 
principle. For instance at the next order the vacuum wave functional reads 
\[/ = e s^J^ J^ where P satisfy the equation 

|^[[r, V], V]l = ^T . P = i I Tr(F,AF,) - e^^'TT{F,[F,, F,]) (122) 

with A the covariant laplacian. Now one can see from the definition of the 

2 

regularised kinetic term that j^T when acting on gauge invariant operators 
of order 6 reproduces a linear combination of such operators plus a sum over 
lower dimensional terms /i^V and /i^ which contributes to a renormalization of 
the vacuum energy and coefficients of V. Thus up to renormalization, which 
can be reabsorbed at this order in a redefinition of M — i> M (1 + -4^), with a 
a numerical constant, we have 

P = j aTT{F,AF,) - be'^'TT{F,[F„ F,]). (123) 

where a,b are numerical coefficients (independent of N,fi,g). These coeffi- 
cients are uniquely computable from the regularization prescription we have 
given. This calculus is quite lengthy but otherwise straightforward and will be 
presented elsewhere. This procedure can be implemented at higher order in a 
strong coupling expansion. 

The conclusion is then that in the infrared regime ( more precisely for 
slowly varying and low amplitude field) and in a strong coupling expansion 
the vacuum wave functional is given by 

^_g-^/Tr(B2)_ (124) 

where the mass scale m computed here only to ffist order in 1/g^N'^ can 
be recursively computed in a strong coupling expansion. That is m = M(l + 
, "j^.2 + ■ ■ ■ ) where the dots refer to terms of higher order in l/{g'^N)'^. We can 
now give easily the physical interpretation of the dynamical scale m entering 
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the definition of tlie infrared vacua wave functional TT. Lets compute the 
expectation value of a large Wilson loop Wc with area Aq 

(Wc) = (^.ac|H^c|^.ac). 

This average should be dominated in the limit Ac — > oo by the expectation 
value with respect to the infrared vacua p24|) and is therefore given by 



{Wc) r^ DAe ^^^^^'^Wc{A). (125) 

Ja/g 

We recognize the partition function of 2+1 Euclidean Yang-Mills theory with 
coupling g2_^_l = ^-^. Now this expectation value reduces in the limit of large 
area to the expectation value of the Wilson line in the 2 + 1 Yang-Mills vacua. 
This vacua has been recently exactly constructed at least in the large N limit 
where no screening due to n-ality is expected fH] and since 2 + 1 Yang-mills 
shows the property of confinement one gets 

{Wc) ~ e^"^^ (126) 

where the string tension is given by 

^ = ^2+1^^- (127) 

Translated back into 3 + 1 Yang-mills coupling one gets 

One should insist here that this results has been only establish firmly only in 
a strong coupling regime at first order, that is m = M, and is not yet shown 
to survive the continuous limit. 

6 Conclusion 

We have seen in this work that is is possible to give a formulation of pure 
QCD in terms of local gauge invariant observables and computed explicitly 
the jacobian for this transformation which shows that the trivial wave func- 
tion is normalisable. We have also constructed the regularised hamiltonian 
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and expressed the vacuum wave function at first order in a strong coupling 
expansion. This expansion can be extended to higher order as an expansion 
in terms of higher dimensional local operators. 

Several problems in this framework have not been explored yet. The most 
formidable is the problem of the continuum limit. As we have seen a mass 
scale makes its appearance in the construction of the vacuum wave functional 
but this scale is cutoff dependent. In order to establish on a robust basis the 
existence of a dynamical mass scale in the infrared one would need to show that 
it is possible to remove the cut-off without destroying the presence of a non 
zero string tension. This would mean first showing that no phase transition 
makes its appearance while going from strong to weak coupling (that is M > 
from strong to weak coupling) and moreover that the scaling of the dynamical 
mass is such that its continuum value is controlled by the usual renormalization 
group scaling [12] 

(47r)2 (47r)2 11 //^^x 34 / //i^XX ^ ^ 

^ ^ ^ ' ' IMV^ +TTln In^ . (129) 



Notice that the detail study of the continuum limit in 2+1 pure QCD from 
the KKN approach needs to be completed too. 

A related project which is less ambitious but phenomenologically more 
interesting concerns the possibility to construct in 3+1 pure Yang-Mills the 
kernel interpolating between the infrared and ultraviolet behavior of the vac- 
cuum wave functional and extract from this an estimate for the mass gap in 
terms of the string tension [Sj. 
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7 appendix 

In this appendix we want to establish the main technical result of this paper, 
that is the proof of the fact that the kinetic term acts diagonally on the 
potential term which is the lowest order local gauge invariant operator (|llip . 
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This amounts, has shown in the main body of the paper, to the computation 
of 

dl,dl{TiAd{^^.{y,z))l=y. (130) 

In order to perform this calculation it is convenient to introduce some conven- 
tions and establish some basic results. First one recall the definition of the 
current Jj 

diHi = JiHi, (131) 

One also introduce an auxiliary current Jj 

diHi+i = Hi+iJi. (132) 

The integrated Bianchi identity HiH2H^ = 1 allows us to express the third 
type of H derivatives in terms of these currents 

diHi-i = —Hi^iJi — JiHi^i. (133) 

Finally, derivatives of the currents J can be purely expressed in terms of J 
derivatives: 

di+iJi = i/j+^SjJj+ii/j+i = H^j^^Bi^iHij^i, (134) 

and 

A-i Ji = Hi_A-iJiH;\ - diJi.i (135) 

The computation of p30|) splits in two parts, first one needs to evaluate 
U{y) ^ [dUdUTT^4HM{dyd!e^{y,z)){HM)-')]^^^ (136) 

dxTr,4[Ci (i/.(l/)(9fGf(y,x))] [Ci id^G^ix,y))Hr\z)]) . 
One starts by the evaluation of 

Ci(^.(2/)(5fGf(y,x)) (137) 

= {dtiHi{y)^i{y, x))6f,{y, x) + Hi{y)Ai{y, x)df_^^6^{y, x) 

where we have used the property (fH^ dfAi{y,x) = 0. Moreover, choosing 
i = 1 without loss of generality and starting from the definitions of Ai{y,x) 
and the currents one can compute 

dlHi{y)Ai{y,x) = Hi{y) 73(2/1,^/2,^3) - ^3(3^1,^/2, 2/3) (138) 

+ H2{xi, 2/2, 2/3) (-^3(3^1, X2, 2/3) - h{.xi,y2, ys)) H2\xi, y2, 2/3)] Ai(2/, x). 
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We can Taylor expand ()137|) and get 

dUH,iy)idfG^,{y,x)) (139) 

= [-X,B,{y) + X^{H^^B,H^){y) + 2^2X3] 5,{X)H,{y) + 0{X^) 

with Xi = Xi — yi and where we have made use of the relations ()132j) . A(x, x) = 
1 and the explicit form of the regulated delta function derivative: 8^5 ^{x, y) = 
2ij,'^{x3-ys)6^{x,y). 

One continues with the evaluation of 

Ci((5fGr(^,l/))^r'(l/)) (140) 

= -Ux,y)Hr\y)dl,6^{x,y) + dl, {dyUx,y))Hr\y)G^{x,y)) 

with 

9f Ai(x, y) = Ai(x, y)H^^{y) [hiyi, 1/2, 1/3) - -hiyi, I/2, x^)] Hi{y) 
+H2{yi,X2,X3) Jiiyi,X2,X3) - Ji{yi,y2,X3) H2^{yi,X2,X3)Ai{x,y), 

and 

dUd!Ai{x,y))H^\y) = A^{x,y)H^\y)d3My). (141) 

We have also used the property 

^|A^{x,y)H^\y) = 0. (142) 

The Taylor expansion of ()140|1 is therefore given by 



dU{H.{yMG^{x,y)) (143) 

= H^\y)[-2fi^X36,{X) + B2iy)Gt{X) 

+ [X3B2 - X2{H^B3H^')]iy)d3G'tiX)] + ©(X^) 

we have denoted G'f(X) = Gl^{x,y) and similarly for 5^. 
Putting these results together one obtains 

h = jdXTi,, {[-X^B2{y) + X2{H,'B,H3){y) + 2/1^X3] 6,{X) 
[-2fi'X36^{X) + B2{y){l - 2/i2x|)G'^(X) + 2fi^X2X3{HiB3H^'){y)Gl{X)]) 

where we have neglected inside the parenthesis in the integrand all the term of 
cubic and higher orders in Xj and express d^Gi = — 2/i^X3Gi. Due to parity 
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symmetry only the terms containing XfS'f^{X), XiGiS^ or X'^XiGiS^ are non 
zero, the corresponding integrals can be explicitly evaluated 

UxXl6l{X) = -f-^, (144) 

J 4(27r)2 

[ dXX,G,iX)6,{X) = 2-^, (145) 

2^? UxxlX,G,{X)5,{X) = -^, (146) 

J 4(27r)2 

2^' UxxlG,{X)G,{-X) = 2-f-^. (147) 

J 4(27r)2 

(the last integral will be used later.) Moreover for the same reason the cubic 
terms which have been neglected do not contribute to the integral, the first 
non trivial contribution starts at quartic order and the corresponding integral 
is of order l//i, which justifies, a posteriori, our approximation. Eventually, 
using the relation TiadiB^) = 2NTt{B^), Tr„rf(l) = A^^ _ ^ ^^^ g^^g 

h = -f-T {-2NTt{BI) - A^,\N^ - 1)) + O {-) (148) 

4(27r)2 V^y 

which proves ()114|) . 

For the second part of the computation one needs to evaluate 



Ih{y) ^ [dUdl{Iv,,[DyDlQl,{y,z))]^^^ (149) 

dxTT,,{[dUD\GU{y,x)][dUD\GU{x,y)]) 

One first focus on the derivative (take i = 2) 

[d\DyG^,{y,x)] (150) 

= {dlDyA,iy,x))Gt{y,x) + iDyA,iy,x))5,{y,x) + A,iy,x)dyS,iy,x), 

with 

DlA,{y,x) = H2{ymH^\y)A,{y,x)) (151) 

= {J2{xi,y2, ys) - J2{yi, 1/2, y3))Ai(y, x) (152) 

and 

d!DyA,iy,x) = -dMy)A,{y,x) = -B^{y)A^{y,x). (153) 

31 



The equality ()151|) follows from the definition of Di and the next one from a 
direct computation. The Taylor expansion of the first term is therefore given 
by 

[dlDlG^,{y,x)] (154) 

= {BM{X^5,{X)-G^{-X)) + 2^?X25,{X)) + 0{X^) (155) 

with Xi = Xi — Hi. 

The second term of p49j) is given by 

[d\DlG^,{x,y)] = {dlDlk,{x,y))G^,{x,y)-{Dlk,{x,y))5,{x,y) 
+ (af Ai(x, y))dlGi{x, y) - Ai(a;, y)d^6^{x, y), 

with 

D2Ai{x,y) = dl{KiH2){x,y)H^\y) (156) 

= H2{yi,X2,x-i)H3{yi,y2,y3)[J2iyi,y2,X3) - J2{y)]Hi{y) 
= X^{H2B^H:^'){y) + 0{X^) (157) 

and 

d\D2Mx,y) = X^[H2{D^B,)H^%y) + 0{X''), (158) 

9iAi(a;,y) = X2BM-X3{Hi'B2H{){y) + 0{X^), (159) 

with Di = di + Ji. Thus, up to terms of order l//x we have 

Ihiy) = jdXTTad{[B,{y){X,6^{X)-G,{-X)) + 2fi^X26^{X)] (160) 

X Xs{H2DiBiH^'Gt{X) - H2B,H^'6^{X) - H,B2Hi2fi^Gt{X)} 
+B, 2fi^XlGtiX) - 2ii''X25,iX)] ) 

Due to parity symmetry the term proportional to X3 do not contribute, we 
are left with 

Ihiy) = TTUBl{y))jdXX',G,{X){X^6,{X)-G,{-X)) 
-Afx^TiUl) J dXXiS.iXy 

= -^ {-2Ni:r{Bl{y)) - A^,\N' - 1)) . 
4(27r)2 

This proves (J114p . 
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